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ABSTRACT 


We  consider  in  this  paper,  pulsating  laminar  flow  superposed 
on  the  steady  motion  of  a  viscous  incompressible  electrically 
conducting  fluid  in  an  annular  channel  between  two  infinitely 
long  circular  cylinders  under  a  radially  impressed  magnetic  field. 
The  general  magnetohydrodynamic  equations  are  simplified  by  the 
conditions  of  the  problem  to  three  equations  in  pressure,  velocity 
and  magnetic  field.  One  equation  gives  the  pressure  variation  in 
the  radial  direction;  the  other  two  are  coupled  equations  for  the 
velocity  and  magnetic  field,  which  are  functions  of  the  radial 
variable  only.  The  solutions  of  these  equations  have  been  obtained 
in  terms  of  Kelvin  and  Lommel  functions  and  the  velocity  profiles  are 
examined  for  different  values  of  the  frequencies  of  pulsation. 


PULSATING  MAGNETOHYDRODYNAMIC  FLOW  SUPERPOSED  ON  THE 
STEADY  LAMINAR  MOTION  OF  AN  INCOMPRESSIBLE  VISCOUS 
FLUID  IN  AN  ANNULAR  CHANNEL 


M.  N.  L.  Narasimhan 


Introduction 

In  recent  years  flows  through  small  pipes  have  been  investigated  with 
increasing  interest  in  the  circulation  system  of  the  blood  and  recently  theories 
of  pulsating  flow  are  applied  to  the  supercharging  system  of  reciprocating 
engines  and  the  surging  phenomena  in  power  plants,  and  so  forth  (Uchida  1956). 

But  it  would  be  of  even  greater  interest  to  an  analyst  or  experimentalist 
to  investigate  similar  types  of  flow  in  an  annular  channel  in  the  presence  of 
magnetic  fields,  since  the  effect  of  impressed  magnetic  fields  on  pulsating 
flows  are  of  considerable  physical  significance  particularly  in  space  technology 
problems  and  also  biophysical  problems. 

Globe  (1959)  has  considered  the  laminar  steady  state  flow  in  an  annular 
channel  in  the  presence  of  a  radial  magnetic  field.  In  the  present  investigation 
we  study  analytically  the  pulsating  flow  superposed  on  a  steady  flow  in  an 
annulus  under  a  radial  magnetic  field.  We  assume  in  the  manner  of  Globe  that 
an  approximation  to  the  desired  radial  magnetic  field  may  be  obtained  by  the 
use  of  a  permeable  core  within  the  annulus  and  a  permeable  cylindrical  shell 
outside  the  annulus.  The  flux  lines  would  close  through  these  permeable  paths 

at  long  distances  from  the  region  of  interest.  The  source  of  the  flux  could  be 
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disks  of  permanently  magnetized  material  between  the  permeable  paths  and  the 
annulus  channel. 

The  present  investigation  aims  mainly  at  the  analysis  of  the  pulsating 
motion  under  such  a  radial  magnetic  field. 

We  now  consider  an  infinitely  long  annular  channel  of  inner  radius  a  and 
outer  radius  b,  as  shown  in  cross  section  in  figure  1. 


Figure  1.  Annular  channel  with 
impressed  radial  field. 

A  radial  magnetic  field  HQ  =  ~  ,  where  a  is  a  constant  is  impressed  across  the 
channel. 

We  shall  now  show  that  the  independence  of  the  electromagnetic  quantities 
with  respect  to  the  coordinate  in  the  direction  of  flow,  usually  assumed  in  plane 
solution,  is  a  necessary  consequence  of  certain  conditions  of  the  problem. 
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Governinq  Equations 

On  the  assumptions  that  (1)  the  fluid  is  incompressible  viscous,  (2)  the 
displacement  current  and  free  charge  density  are  negligible,  (3)  the  permeability 
and  conductivity  are  constant  scalar  quantities,  and  (4)  the  Lorentz  force  is 
the  only  body  force  on  the  fluid,  the  magnetohydrodynamic  equations  (Cowling 
1957)  are 


VX  H  =  J  , 


p  iLM. 

VXE  =  , 


V-  H  =  0 


J  =  o-  (E  +  pV  X  H  )  , 


v-v  =  0  , 


3V  -*•  — ► 

g  +  (V.V)V 


-Vp  +  pvV2V  +|JiTxH  . 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


In  these  equations  H  is  the  magnetic  field,  E  the  electric  field,  p  the 
permeability,  J  the  current  density,  V  the  velocity,  p  the  pressure  and 
v  the  Kinematic  viscosity.  (M.K.  S.  units  are  understood). 

We  obtain  now,  in  addition  to  equation  (6),  one  more  coupled  equation  for 
H  and  V  by  eliminating  J  and  E  among  (1),  (2)  and  (4)  .  Making  use  of 
the  divergence  relations  (3)  and  (5),  one  obtains  the  following  equation: 

-  (H-V)  V  +  (V -V)  H  =  —  V2  H  . 


(7) 
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Equations  (6)  and  (7)  must  now  be  expressed  in  cylindrical  coordinates 

r,  0,  z  .  The  equations  (6)  and  (7)  simplify  greatly  because  of  the  following 

conditions: 

0 

(1)  We  assume  axial  symmetry  around  the  axis  so  that  —  =  0  . 

(2)  Since  the  flow  is  laminar  and  parallel  to  the  axis,  v  =  v  =  0  . 

9v 

2 

Furthermore,  from  ( ^ ,  ~r—  =  0  . 

O  Z 

(3)  We  assume  that  the  applied  field  HQ  =~  fixes  the  normal  component 

of  the  magnetic  field  at  r  =  a  and  r  =  b  for  all  values  of  z,  and  that  this  is  the 
only  field  impressed.  Several  consequences  follow  from  this  assumption: 

(a)  from  the  z-component  of  (4), 


=  a 


Ez  +  jx(VXH}z 


(8) 


we  note  that  j  must  vanish.  This  is  because  E  can  arise  only  from  an 
z  z 

applied  E  field  or  free  charges  in  the  flow,  neither  of  which  exists,  and  because 
(V  XH)z  must  vanish,  since  V  has  a  z-component  only.  Similarly  from  the 
r-component  of  (4)  j  vanishes,  (Globe  1959).  Since  H.  can  arise  only  from 
an  r  or  z  component  of  current,  there  can  be  no  0  component  of  H  from  the 
currents  in  the  fluid.  Since  there  is  also  no  impressed  magnetic  field  in  the 
6-direction,  it  follows  that  H  vanishes  everywhere  in  the  channel. 

U 

(b)  In  the  manner  of  Globe  (1959)  let  us  consider  a  section  of  channel 
bounded  by  transverse  planes  at  z  =  ±h  (figure  2).  Applying  the  divergence 
theorem,  and  Equation  (3),  to  this  volume,  we  get 
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Iff  (V  •  H )  dQ  =  //  H  •  ds  =0  . 
Volume  Surface 


But 


//  H  •  ds  =  -2  7T  ah  Hf(a)+  2:rb  Hr(b) 

+  [TTb2  -  TT-a2]  [HzCh)  -  Hz(-h)]  .  (10) 


A 

777777777777777 

b 

iLa _ 

- -  - - - 

mjTTinzjznjjJi 

Z  =  -h  Z  =  h 

Figure  2.  Section  of  channel  bounded  by 
transverse  planes. 


Since  the  first  two  terms  on  the  right-hand  side  of  (10)  cancel  each  other,  it 

follows  that  Hz(h)  =  Hz(-h)  .  Since  h  is  arbitrary,  we  obtain  (9Hz/0z)  =  0 

This  conclusion  is  independent  of  the  nature  of  flow. 

9H 

(c)  Since  =  0,  it  follows  from  (3)  that  -^-(riy  =  o  .  The 
radial  field  must  then  be  equal  to  the  impressed  field,  and  is  unaffected  by  the 
flow  and  independent  of  z  . 
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With  these  simplifications,  the  equations  (6)  and  (7)  yield: 


a  8H 

If  +  l*H,  z 


8r 


z  8r 


0  , 


(11) 


9v  8H 

2  _ii£_z  =  _££  +  v 

p  8t  r  8r  9z  p 


1  82v 


\  8r 


,  8v 

_z  1_ _ z 

2  +  r  8r 


(12) 


8H  9v  8  H 

z  uatr  _ z  _ 


.  8H 

_ _ z  1 _ z_ 

^  8t  "  r  8r  =  .  2  +  r  8r 

8r 


(13) 


8p 

H  and  v  are  functions  of  r  and  t  only.  It  follows  from  (12)  that  — - 
z  z  oz 

must  be  independent  of  z  .  By  differentiating  (11)  with  respect  to  z,  it  can 

be  seen  that  is  independent  of  r  also.  Hence  -~f  must  be  a  function  of 

time  t  only.  We  may  therefore  set 


-9J?=  _ 

9z 


P(t) 


7 


where  P(t)  is  a  function  of  t  only. 

Once  Hz  is  determined,  the  variation  of  p  across  the  channel  may  be 
found  by  integrating  (11): 


pHz 

p(r,  z,t)  +-J~  =  -  P(t)z 


Now  (12)  and  (13)  can  be  rewritten  as: 

=  P(t)  +  pv 


9v  9H 

_ z  jia  z 

9t  "  r 


9r 


9r 


(14) 
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fA(T 


3H  9v 

_ z  uacr  _ z 


at 


3r 


a2tt  ,  8H 

_ z  I _ z 

Q  2  r  ar 

8r 


(15) 


The  boundary  conditions  for  our  problem  are: 

Vz(a,t)  =  0,  Vz(b,t)  =  0  ,  .  (16) 


H  (b,t)  =  0  ,  (17) 

z 

3H 

*5p(b,t)  =  0  .  (18) 

Equations  (16)  contain  the  no  slip  conditions.  Equation  (17)  follows  from  the 
fact  that  J  has  a  0  component  only,  so  that  the  currents  in  the  annular 
channel  are  like  those  in  an  infinite  solenoid.  These  currents  will  therefore 
produce  no  field  for  r>  b,  and  since  there  is  no  impressed  field  in  the  z-direction, 
continuity  of  the  tangential  component  of  H  requires  (17)  to  be  true.  Equation 
(18)  is  obtained  as  follows.  We  have 

3H 

.  _  _ _ z 

J  0=  "  9r 


But 


and  V  must  vanish  at  r 


o-(i(VXH)e  , 

Hence  j  must  also  vanish  there  and 
0 


3H 
_ z 

3r 
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The  Solution 


Let  us  introduce  the  following  non-dimensional  quantities  and  parameters 


X  ,  T  = 


tv 

2  -  9 
a  (p.<rv)2 


(v/a)  ’  (v/a)(o-Pv)2 


.  3  -1  -Z  do  3-1  -2 

$(t)  =  -apv  f-  =  a  p  v  P(t)  , 


,  2  2 

u.  a  cr  ,  .? 

P  =  ^7“  »  v  =  (n°-v)2  , 


and  let  —  =  5  ; 
a 


p  is  a  form  of  Hartmann  number  appropriate  to  an  annular  channel.  On 
introducing  these  non-dimensional  variables  and  parameters  into  the  equations 
(14)  to-tlS),  we  obtain  the  non-dimensional  form  of  the  latter  system; 


JL  9V  _  9_V  1  . 

y  9t  “  8V2+  x  ax  +  \  ax  +  ^T) 


9H  _  9_H  1  9H  £  9V 

9t  „  2  +  x  ax  +  x  ax 


V(X  ,  t)  =  0  at  X  =  1 
V(X  ,  t)  =  0  at  X  =  6 
H(X  ,r)  =  0  at  \  =  6 


0H 

0^(X  ,t)  =  0  at  X  =  6  . 
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The  non-dimensional  pressure  gradient  4>(t)  is  a  function  of  non-dimensional 
time  t  only  and  hence  we  can  express  it  in  the  following  form 

4>(T)  =  +  ^1  ®  (25) 

where  and  4>^  are  constants  and  may  be  assumed  to  be  real  for  the  sake  of 
simplicity  and  u  also  a  real  number  denoting  the  frequency  of  the  vibration. 

The  corresponding  solutions  for  V  and  H  will  now  be  assumed  as: 


V(^,T)  =  V0(V)  +  V^M  eiWT  (26) 

H(V,t)  =  HQ(X)  +  HjCM  eiwT  .  (27) 


The  solutions  for  <j>,  V  and  H  are  respectively  the  real  parts  of  the  expressions 
in  (25),  (26)  and  (27). 

Substituting  these  forms  (2$,  to  (27)  into  equations  (19)  to  (24)  we  have  the 
following  sets  of  equations: 


ifo 

,  2  +  X  dX  +  X  dX  +  *0  °  » 

aX 


(28) 


^0  1  2,  +  £?i 

dx  2  +  X  dX  X  dX 


=  0 


(29) 


I 
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and  VQ  =  0  at  X  =1,6 
HQ  =  0  at  X  =6 


a 

dX 


0  at  X  =  6 


(30) 


.  d2V.  ,  dV,  dH,  . 

and  - l+L_l+£__l  +  4)  ,  M.  v 

2  r  X  dX  X  dX  91  y  V1 


(31) 


..  2  +  X  dX  +  X  dX 
dX 


and  Vj  =  0  at  X  =  1, 6 


Hj  =  0  at  X  =  6 
dH 


*  IwyHj  , 


(32) 


_ 1 

dX 


=  0  at  X  =  6 


(33) 


First,  we  shall  determine  the  steady  part  of  the  solution  from  the  set  (28)  to 
(30). 

dH0 

We  integrate  (29)  once  and  substitute  for  into  (28);  then  we  solve  the 

resulting  second  order  equation  for  VQ;  and  then  we  substitute  the  resulting 
expression  for  VQ  into  the  first  integral  of  (29),  and  finally  integrate  once  again 
to  get  Hq.  We  obtain: 
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sinh  (P  In  6} 


J  P  *4  ,  (34) 


^  4 

V„(\)  =  ? —  \2  In  X  -64X2ln^-  -  In  6  ;  p2  =  4  , 

4(6  -1)  I  °  X 


r  2  X, 

-p*Q  ^2  &2  6  (icosh(pinX) -cosh(pin6)}  -cosh(pin— )+ 1 

H0^  =  aZ  A  2  “  p  sinh  (p  In  6) 

p  -4  L 

P2*4  (3 


..  .  *0  .4.  2  X.  64  .  X. 2  lnX.  .2.  .  (64  -  1V(62  -  X.  2) 

’  H0(X}=“T7;  6X  lnT-  —  ln6— Y~f6  In 6  +  * 

2(6  -1)  L  2X. 

P2  =  4  (37) 


Next  we  consider  the  equations  (31),  (32)  and  (33).  The  solution  of  this 

i 

system  of  equations  for  general  values  of  y  =  (jicrv)2  runs  into  difficulties. 
Hence  to  deal  with  the  actual  physical  situations,  it  will  be  sufficient  to  solve 
the  system  for  y<<  1,  since  for  most  of  the  incompressible,  electrically 
conducting  viscous  fluids  on  the  surface  of  the  earth,  y  <<  1  . 

For  instance: 


Fluid 

Hg(20°c) 

Na(500°c) 


1.55  X  10 


0.47  X  10“3 


Pb(500°c ) 
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Hence  we  shall  assume  y  k<  1  and  ~  y  •  Then  writing  k  =  JT  ,  we 

2 

find  that  the  equations  (31),  (32)  and  (33)  after  neglecting  terms  of  0(y  ),  reduce 
to: 


d25tidA+ifi,0 

J(2f\  d\  1  di 
d  x 


(38J 


d2v 


dV 


dH 


u  ¥1  1  V1  fl  1  2 

- -1  +  —  — -  +  ^  — -  +  A  -  Dc  V  =  0 

J%  2  +  X  dX  +  X  dX  +*1  **  V1 
d\ 


(39) 


with  the  conditions 


V,  =  0  at  X  =  1,  6 


H!  =  0 


at  X  a  6 


(40) 


dHl 

“dX 


=  0 


From  (38)  after  integrating  once  and  using  (40)  we  obtain  • 


a. -A  . 

dX  p  X 


(41) 


Substituting  (41)  into  (39)  we  obtain 

dV 


dV 


2  1  1  ?  2  7  7 

K  - 2  +  X  "dX  "  (ik  X  +  P  }  V1  =  "X\  > 

dX 


(42) 


which  is  a  non-homogeneous  modified  Bessel  equation  with  complex  argument. 


d 
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The  general  solution  of  (42}  is  therefore  of  the  form 


Vj(M  =  Aj  Ip(kXi*)  +  A2  Kp(kXi*)+~  Sj  ^kXi^}  , 


(43) 


where  functions  Ip,  Kp  are  modified  Bessel  functions  of  the  first  and  second 
kind  respectively  and  of  order  p  ,  and  ^  is  a  Lommel  function  (Watson 
1944)  which  occurs  as  a  particular  integral  of  (42) .  A^  and  are  arbitrary 
constants.  This  particular  integral  can  be  written  either  as  ascending  or 
descending  power  series  in  the  argument  (Watson  1944)  according  as 
1  *  P  *  “  (2p-  1)  or  1  *  P  = -(2p- 1),  where  p  is  a  positive  integer.  For  the 
sake  of  simplicity  we  shall  demonstrate  the  solution  when  1  *  p#  -(2P-1)  . 
This  condition  is  satisfied  when  p  =  an  odd  positive  integer  .  Thus  we  shall 
demonstrate  the  solution  when  P,  the  Hartmann  number  is  an  odd  positive 
integer  say  equal  to  2q  -  1,  q  being  positive  integer. 

Hence  in  this  case  we  obtain  the  complete  solution  of  (38)  as: 


{K3(ki*)IB(k\iVla(ki‘)  Ka(k\i^)> 


*1 


'S1  {Kp(k6i^Ip(k\i?)-Ip(k6i^)Kp(k\i2}} 

— ? -  - +  s^kXi2) 


IJki2)  KJk6i2)  -  IJk6i2)  K  (ki2) 


(44) 
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where  the  <XR#;reeslon  for  the  Lommel  function  Is 


.  .  t  z 

Si  S  \l)  -  _  J  -  7  ■>  77 

^  *  (m+1)  -v  {(|a+1)  -v2  }  {(,*+  3)2-v2> 


and  p.  *  vj"'  -t^ap-O,  P  being  positive  integer. 

The  fuiooMtioans  and  Kp  can  also  be  expressed  in  terms  of  Kelvin  functions 
or  in  termsoflocfth_eir  omoduli  and  phases.  For  real  x,  we  have 


Ip(xi  5| »  =e  '^(beyc  +  ibeipx)  =  Mp(x)  exp  {i[0p(x)  ~~  M) 
KpCxiV  =e^  (kerpx  +  ikei^x)  =  Np(x)  exp{i[i|ip(x)  +  j  p*]  } 


where  Mpai6and  Np  are  their  respective  moduli  and  0p  and  i}jp  are  their 
respectivepHqphas  es.  In  an  analogous  manner,  we  shall  write 

n  i. 

Sl,  '  '  =  ppCx)4-  iQp(x)  =  Lp(x)  exp  {i[Xp(x)]> 
where  Lfi(i|  (*HS  the*  modulus  and  XpM  is  the  phase  of  s.  (xi^"},  P  and  Q 

r  r  *>P  P  P 

being  the  leaesal  a.  nd  imaginary  parts  of  s  (xi^J. 

A  >  r 

Once  theffhe  s  olutdon  for  Vj(\)  is  known,  we  can  now  determine  Hj(X)  from 
(41).  We  lavMve  *hus 
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iP-t*! 


’ pCkei^)  {Kp(ki^)  [f (x )  -  f(6)]  - yki*)  [g  (x }  -  g(6)] } 

-s1)p(kl^){KQ(kSi^[f(X.)-f(6)3-IQCk61^)[g(X>-g(6)] 
IpCki^J  KpCkSi1 )  -  Ip(k6i*  }  Kp(ki*  ) 


+  h(X }  -  h(6) 


(45} 


where 


X  I  (kXl2) 

f  (X }  -  J  - dx  , 

\  K  (kXi2) 

g(X)=  /  -  dX  , 


X 

h(X)  =  / 


Sl,pCkXi^ ) 

X 


dX.  . 


Thus  the  velocity  distribution  and  the  magnetic  field  are  now  given  by  the  real 
parts  of  (26)  and  (27)  where  VQ(X},  V^X),  HQ(X)  and  H^X)  are  given  by 
(34}  to  (37)  and  (44)  to  (45). 

We  next  obtain  the  real  parts  of  those  expressions  after  expressing  the 
modified  Bessel  functions  and  Lommel  functions  in  terms  of  their  moduli  and 
phases.  Thus  if  we  write  the  velocity  distribution  as 

V(X,T)  =  V0(X)+ Vt(X,T)  , 


where  Vt(X,r)  denotes  the  time-dependent  part,  then 


V,(X,T|=- 


x  y  -x  y 

- p-+  L  (kX)  sin[xa(kX)J  ^  cos  o>t 

l  +  y1 


} 


X1X2  +  yly2 

+  <  4,j.  *  £  +  L  (kX)  cos  [XR(kX)I  )  sin  wr 

xj  +  yf  p  p 


(46) 
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where 


V 


G 


k,k6  k6,k 


C°Sak6,k  C0Sak,k6 


V 


ak,k6  k6,k 


Sinak6,k  Sinak,k6 


X2 


Y2  = 


GkX  ,k,k6 

Gk,kX,k6 

GkX  ,k6,k 

Gk6,kX,k 

GoS  ak,kX  ,k6 

CosakX,k,k6 

CoS  ttk6  ,kX  ,k 

CoS akX,k6,k 

GkX  ,k,k6 

Gk,kX ,k6 

GkX ,k6,k 

Gk6,kX,k 

Sinak,kX,k6 

Sin°kX,k,k6 

Sin“k6,kX,k 

Sin°kX,k6,k 

Gx,y  *  VX)  VY) 

*x,y  =  VX)+¥Y) 
Gx,y,z  =  Mp(x)  Np(y)  Lp(z) 
and  ax,y,  z  =  0p(x)  4/p(y)  xp(z)  • 
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j  2  2 —  1  beinX 

M  (x)  »  vber  x  bei  x  ;  0  (x)  =  tan’1  t—* 


bV 


I  2  2  1  keiaX 

Np(x)  =  v  kerp  x  +  keip  x  ;  ^(x)  =  tan  ’ 


kv 


J~2 - ? -  ,  Qa(x) 

Lp(x)  =  s  Pp  (x)  +  Qp  (x)  ;  Xplx)  =  tan  ' 


and  P  (x)  ,  Qp(x)  being  respectively  the  real  and  imaginary  parts  of 
have  the  following  expressions 


i 


(xi2  ) 


w  = 


(22  -  p2}  (42  -  p2)  (22  -  p2)  (42  -  p2)(62  -  p2)  (82  -  p2) 


4*  •  •  •  • 


OM  =  -  ,  ,  t— — — - — - — _ 

P  22-p2  (22  -  p2)(42 -P2H62  -  P2) 


Total  mean  volumetric  flow  G,  flowing  in  the  z-direction  is  given  by 


2ir  6 


G  =  i 7  /  dr  /  2rr  \  V(K,t)  dX. 


(47) 


or  G  = 


2H, 


t- 


6-1  14  4  2 

4+2  (2(6  “1}~P(6  +l)coth(pin6}  +  2p6  cosech(p  In 6)} 

P  “  4 


and  '  G  = 


2M>r 


4(6-1)  L 


(64-ll2  4..  .2 

]£,  ”  6  (ln  6) 


when  p  #  4  ; 


when  p  =  4  . 


r 
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Discussion  of  the  Solution 


We  now  Investigate  the  nature  of  the  velocity  distribution  Vt(K  ,  t)  and 

the  influence  of  the  magnetic  field  on  it.  Figures  3  and  4  give  the  time-dependent 
Vt(V  ,r) 

velocity  profiles  — r -  in  both  non-magnetic  and  magnetohydrodynamic. 

*0 

flows.  To  illustrate  the  solution  we  have  chosen  co  =  y  ,  since  co  ~y  >  i. e. , 
k  =  1,  radius  ratio  of  the  annulus  —  =  2  and  «(>.=•  1<|>A  .  The  velocity 
profiles  are  obtained  for  different  values  of  cot,  viz  0°,  to  360°  .  In 
Figure  3,  the  velocity  profiles  are  obtained  for  the  Hartmann  number  p  =  0  and 
cor  =  G*  to  360*  .  In  Figure  4  the  velocity  profiles  are  obtained  for  the 
Hartmann  number  p  =  5  and  cot  =  0°  to  270*  .  On  comparison  of  Figure  3  with 
Figure  4,  it  is  found  that  the  flattening  of  the  profile  in  Figure  4  which  is 
characteristic  of  magnetohydrodynamic  flow  under  a  transverse  field  is  evident. 
Furthermore,  the  magnetohydrodynamic  flow  profiles  get  even  flatter  and  flatter 
as  cot  increases  from  0*  to  90*  .  Then  with  further  increase  in  cot,  the 
flow  is  reversed  in  direction  due  to  pulsatory  motion  with  the  profiles  becoming 
less  and  less  flat  between  cot  b  90*  and  180*  .  With  further  increase  in  cot 
from  180*  to  270*  ,  once  again  the  profiles  get  flatter  and  flatter.  When  cot 
is  between  270*  and  360*  ,  the  original  direction  of  flow  is  restored  back  and 
after  this  full  period  of  360*  for  cot,  once  again  the  profiles  get  flatter  and 
flatter  as  cot  begins  the  next  cycle  of  values. 

A  more  general  problem  of  unsteady  magnetohydrodynamic  flow  in  an  annular 
channel  is  presented  in  another  paper. 


Yt(K*) 


Figure  3 


v,  (x,^) 
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